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ABSTRACT 

Numerical simulations of the stochastic end stage of planet formation typically begin with a pop- 
ulation of embryos and planetesimals that grow into planets by merging. We analyzed the impact 
parameters of collisions leading to the growth of terrestrial planets from recent A-body simulations 
that assumed perfect merging and calculated more realistic outcomes using a new analytic collision 
physics model. We find that collision outcomes are diverse and span all possible regimes: hit-and-run, 
merging, partial accretion, partial erosion, and catastrophic disruption. The primary outcomes of 
giant impacts between planetary embryos are approximately evenly split between partial accretion, 
graze-and-merge, and hit-and-run events. To explore the cumulative effects of more realistic collision 
outcomes, we modeled the growth of individual planets with a Monte Carlo technique using the dis- 
tribution of impact parameters from A-body simulations. We find that fewer planets reached masses 
> 0.7-MEarth using the collision physics model compared to simulations that assumed every collision 
results in perfect merging. For final planets with masses > 0.7ME ar th, 60% are enriched in their core- 
to-mantle mass fraction by > 10% compared to the initial embryo composition. Fragmentation during 
planet formation produces significant debris (~ 15% of the final mass) and occurs primarily by erosion 
of the smaller body in partial accretion and hit-and-run events. In partial accretion events, the target 
body grows by preferentially accreting the iron core of the projectile and the escaping fragments are 
derived primarily from the silicate mantles of both bodies. Thus, the bulk composition of a planet 
can evolve via stochastic giant impacts. 

Subject headings: accretion — Solar System: formation — planets and satellites: formation — Earth 



1. INTRODUCTION 

The end stage of terrestrial planet formation is defined 
by the final competition between planetary embryos as 
they grow via collisio ns into a final set of stable planets 
(e.g., [Chambers 2004|). The last giant impact is invoked 
to explain major differences between the solid planets in 
our Solar System. The large core mass fraction for Mer- 
cury could be the result of a single mantle-stripping im- 
pact event ()Benz et al.| [2007). A planet's spin orientation 
and period is st rongly influenced by the last giant impact 
(|Lissaueri ll993). Giant collisions are also invoked to ex- 
plain the formation of E arth's moon (Canup fc AsphaueJ 
1200 ID . t he Pluto system (ICanup|l2005l) . and the Haumea 
system (jLeinhardt et al.ll2010( l. 

The growth of planetary embryos into planets is com- 
monly modeled using A-body techniques in order to cap- 
ture the details of the multiple gravitational interactions. 
Because typical encounter velocities between embryos are 
in the range of 1 to 3 mutual escape velocities (e.g., 
lAgnor et al.lll999t ). almost all A-body simulations have 
assumed that the two colliding embryos merge into a sin- 
gle body with the same total mass. However, hydrocode 
simulations of oblique collisions between embryos found 
that in elastic collisions, k nown as hit-and-run (HfcR ) 
events (|Agnor fc Asphaud 12001 lAsphaug et"aH l2006h . 
are a common outcome at speeds just above the mu- 



tual escape velocity. Asphaug (|2010f) estimated that 
about half of giant impacts are hit-and-run events. The 
physics of giant impacts is sufficiently complicated that 
the impact parameters demarcating the transition be- 
tween merging and hit-and-run outcomes was defined 
only recently by an empirical an alysis of over a thousand 
hydrocode impact simulations ( Kokubo fc Gendal 120101 ; 
IGenda et al1l2012lh 

Collision outcomes other than merging and hit-and- 
run are possible during planet formation. The range of 
outcomes include partial accretion, partial erosion, and 
catastrophically disruptive events. Some planet forma- 
tion simulations, focusing on specific stages of planet 
growth, have directly modeled each encounter between 
two bodies i n order to capture the details of the colli- 
sion physics (IGenda et al.ll2011l:lLeinhardt &; Richardson! 
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2005^ Leinhardt et al.l 120091 ) . but this approach is too 
computationally expensive to be widely adopted. To 
date, planet formation studies have not had access to a 
robust collision physics model that could be implemented 

easily into A-body simulati ons. 

In a companion paper (jLeinhardt fc Stewart! I2012h . 
we have developed a comprehensive analytic model to 
predict the outcome of any collision between gravity- 
dominated bodies. The model predicts the transitions 
between collision regimes and the size and velocity dis- 
tribution of fragments. Collision outcomes depend on 
the projectile-to-target mass ratio, impact angle, impact 
velocity, and two material parameters that define the 
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catastrophic disruption criteria. Catastrophic disruption 
is defined by the specific impact energy needed to dis- 
perse half the total mass. For collisions that result in 
partial accretion or erosion of differentiated bodies, we 
estimate the chan ge in the core - mantle mas s ratio fol- 
lowing the work of iMarcus et all (|2009l 120101 ). The new 
analytic model now allows planet formation simulations 
to capture the full diversity of collision outcomes. 

To assess the effects of including more realistic colli- 
sion physics in planet formation studies, we analyzed the 
impact parameters from recent iV-body simulations of 
terrestrial planet formation that assumed perfect merg- 
ing and calculated the more realistic collisio n outcomes 
(jO'Brien et al.l [200l iRavmond et~aH [20091 ). The re- 
sults provide insights into the most important collision 
regimes, quantify the probabilities of different outcomes, 
and identify future challenges in planet formation simu- 
lations. 

Then, using the distribution of impact parameters from 
the iV-body simulations, we modeled individual planet 
growth with a Monte Carlo technique in order to assess 
the cumulative effects of diverse collision outcomes. We 
compare the distribution of final planet masses from sim- 
ulations with different assumptions about collision out- 
comes and tracked changes in the core mass fraction. 
The collision physics model allows for more detailed in- 
vestigation of the evolution of the bulk chemistry and 
for tighter isotopic constraints on terrestrial planet for- 
mation studies. 

2. THE COLLISION OUTCOME MODEL 

The outcome of a collision in the gravity regime was 
predi cted using the new analytic collision physics model 
from iLeinhardt fe Stewart! (|2012f ). hereafter LS12. The 
model is summarized here; refer to LS12 for the full de- 
tails and limitations. In the gravity regime, erosive out- 
comes are primarily controlled by overcoming gravita- 
tional forces rather than overcoming material strength; 
bodies larger than ab out 100 to 1000 m in s ize are domi- 
nated by gravity (e.g..lMelosh fc Ryan 1997|). The model 
predicts the collision outcome as a function of projectile- 
to-target mass ratio (7 = M p /M t ), impact angle (6), 
and impact velocity (V*). By convention, the projectile 
is less massive than the target (M p < M t ). The follow- 
ing equations were used to generate the boundaries and 
contours on the example collision outcome maps shown 
in Figure [IB- 
Partial accretion and erosive outcomes are defined by 
the catastrophic disruption criteria (Qrjj), which is the 
specific impact energy required to disrupt and gravita- 
tionally disperse half the total mass. For a particular 
collision, the total specific impact energy in the center of 
mass frame is given by 



(fl99l : 



Qr = 



2M tc 



where M tot = M p + M t and /i = M p M t /M tot is the 
reduced mass. For bodies in the gravity regime, the head- 
on disruption criteria increases wit h increasing mass and 
follows the general form derived by iHousen fc Holsappld 

1 A computer code to calculate individual collision outcomes and 
generate collision outcome maps is available from the authors. 
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where Rci is the spherical radius of the combined projec- 
tile and target masses at a density of p\ = 1000 kg m~ 3 
and G is the gravitational constant. At the critical im- 
pact velocity (Vi — V*), the specific impact energy is 
equal to the catastrophic disruption threshold (Qr = 
Q*rd) f° r the particular impact scenario (Rci and mass 
ratio). Instea d of target radius, th e varia ble Rci was 
introduced by Stewart fc Leinhardt] (|2009f ) in order to 
generalize the dependence of Q* RD on size for different 
density bodies and for different projectile-to-target mass 
ratios. The coefficient q g is of order 1 and defined below. 
p, is the velocit y exponent in the point sou rce coupling 
parameter from IHousen fc H olsapplc (|1990l ). LS12 fit a 
value of j2 = 0.36 ± 0.01 using the results from several 
numerical studies of disruption in the gravity regime on 
a wide variety of materials (weak rock, ice, rubble piles). 

LSI 2 derived the disruption criteria for a reference 
case of a head-on impact between two equal-mass bod- 
ies, which represents the minimum energy scenario to 
disrupt a particular target. The critical impact velocity 
for head-on equal mass collisions is given by 



V* 



32nc* 



1/2 



(p^y^Rd 



(3) 



where c* is a nondimensional material parameter that 
is fitted to simulation data. LS12 fit c* = 1.9 ± 0.3 for 
fluid planets and c* = 5 ± 2 for various solid planetesimal 
analogs using the results from several numerical studies, 
c* is defined such that it represents the value of the equal- 
mass head-on disruption criteria in units of the specific 
gravitational binding energy of the combined mass: 



Q 



RD tl =l 



c*-irp 1 GR cl . 
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(4) 



Equation [5] is named the "principal disruption curve." 
Combining the previous equations gives the coefficient 
for Equation [21 



1 /32nc* 



1 9 = ~ 8 



3/3/2 



(5) 



The derivation of Equation [2] assumed that all of the 
projectile's kinetic energy is deposited in the target. This 
assumption is often invalid as a portion of the projectile 
misses the target for many combinations of projectile- 
to-target mass ratio and impact angle. For an oblique 
impact, LS12 derived the mass fraction of the projectile 
that geometrically overlaps with the target, 

if R t > b{R p + R t ) + R p , 
[p(ttR p 1 2 — 7rZ 3 / 3) otherwise, 

where p is the density of the projectile, b = sin8 is the im- 
pact parameter (6 = for head-on impacts) and l/(2R p ) 
is the fraction of the projectile diameter that overlaps 
the target. The impact angle and overlapping mass frac- 
tion are defined at the point of first contact assuming 
spherical bodies. Using just the interacting fraction the 




Figure 1. Predicted collision outcome maps using the analytic model for strcngthless planets (with material parameters c* = 1.9 
and p. = 0.36) for two projectile-to-target mass ratios of (A) M p : Mt = 1 : 10 and (B) 1:40. Impact velocity is normalized by the 
mutual surface escape velocity; impact parameter is spaced according to probability. Colored regions denote perfect merging (dark 
blue), partial accretion to the target (medium blue), net erosion to the target (white), hit-and-run (green), and graze-and- merge 
(blue-green). Vertical red line at b cr ; t denotes the transition between non-grazing and grazing events. Thick black curve - critical 
velocity for catastrophic disruption; dashed grey curves - 10% and 90% of target mass in largest remnant; dotted black curve - 50% of 
projectile accreted; dot-dashed blue curve - catastrophic disruption of the projectile in a hit-and-run; das hed blue curve — o nset of erosio n 
of projectile in a hit-and-run. Symbols denote proposed giant impact events: Mercury ( Bcnz ct al. 2007); ■- Earth-Moon ( Canup 2004). 



projectile, the angle-adjusted reduced mass is 



Mo 



aM p M t 



(7) 



The disruption criteria and critical impact velocity for 
a specific mass ratio and impact angle are calculated by 
adjustments to the head-on equal-mass principal disrup- 
tion curve. To achieve disruption with a smaller projec- 
tile, the impact velocity rises to 

N2T V(3a0 



v* = 



1(7+1)' 
4 
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V* 1 



(8) 



and, from Equation |2J the catastrophic disruption crite- 
ria increases by 



2-3/3 



Qrd — Q 



RD,j=l 



'AD,7=1 




2/(3/i)-l 



(9) 



Next, accounting for the projectile interacting mass frac- 
tion for oblique impacts, 



Qrd — ~ 



2-3/1/2 



Ql 



RDi 



(10) 



where the prime notation indicates the value for an 
oblique impact. Finally, from Equation [I] the critical 
impact velocity adjusted for both the mass ratio and im- 
pact angle is 



V* = J2Q'* D M tot /fi 



(11) 



The equations for Qr D and V* are curves as a function 
of Rci, impact angle (within fi a ), mass ratio (within /i), 
and two material parameters, c* and /2. 



The specific collision outcome regime is determined us- 
ing the mutual escape velocity, disruption criteria, and 
impact parameter. When the impact velocity is be- 
low the mutual escape velocity of the interacting mass 
(M t ' ot = aM p + Aft), the two bodies are assumed to 
merge completely (perfect merging). Above V^,' sc = 
yj2GM[ ot j (Rp + R t ), the collisions outcomes are divided 
into two groups: grazing and non-grazing. 

The transition between non-grazing and grazing out- 
comes is d emarcated by a critical impact parameter, & C rit- 
Following Asphaug (20lJJ), the center of the projectile is 
tangent to the target at the critical impact parameter, 



b c ,r\t 



R 



R + r 



(12) 



Non-grazing (b < 6 cr ;t) collisions transition from perfect 
merging to the disruption regime with increasing impact 
velocity. In the disruption regime, the outcome may be 
partial accretion or partial erosion depending on the mass 
of the largest remnant, Af| r . The largest remnant mass is 
calculated using the impact energy and the catastrophic 
disruption criteria. For collisions with < Qr/Q% d < 
1.8, Mi r is proportional to impact energy: 



Aflr 

Aftot. 



= -0.5 



Qr 

Q 



- 1 



RD 



0.5. 



(13) 



Because a single line fit a wide variety of simula- 
tion results, we named this linear relationship the 
"universal law for the mas s of the largest remnant" 
(| Stewart fc Leinhardti 120091) ■ In the case of super- 
catastrophic disruption, Qr/Q'rd > 1.8, the largest rem- 
nant follows a power law, 



Afi r 

Mint 



0.1 

L8^ 



Qr 
Qrd 



(14) 



where r) ~ —1.5 based on laboratory disruption experi- 
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ments, as discussed in LS12. 

Grazing (b > 6 C rit) collisions transition from perfect 
merging to a graze-and-merge regime, where the two 
bodies hit and separate (cleanly or not ) but are gravi- 
tationally bound and ultimately merge ([Leinhardt et all 
I2010T) . The boundary between graze-and-merge and hit- 
and-run, where the two bodies hit and then escape each 
other with the target relatively intact, is difficult to de- 
fine. From over a thous a nd hy drocode collision simula- 
tions, [Kokubo & GcndaJ (|2010| ) fit an empirical relation 
for the velocity of transition to hit-and-run, Vh r , 



= Cl ^(l ^ bf" + c 2 ( z + c 3 (l - bf' z + c 4 , (15) 

t'csc 

where ( = (M t - M p )/M tot , c x = 2.43, c 2 = -0-040 8, 
c 3 = 1.86, and c 4 = 1.08 (see also iGenda et al.l (j2012ft l. 
We stress that the transition between graze-and-merge 
and hit-and-run is likely to depend on material proper- 
ties of the planet. These coefficients are fit to smoothed 
particle hydrodynamics simulations of fluid Earth-mass 
planets. 

In iV-body simulations, bodies are described by their 
mass, and radii are usually estimated assuming a 
constant bulk density in order to calculate the im- 
pact angle and mutual surface escape velocity, V^ sc = 
y/2GMtot/{R v + Rt)- A bulk density of 3 g cm~ 3 was 
used in the two studies considered here, so we adopted 
the same value to calculate R p and Rt in this work. 

For grazing collisions with impact velocities exceeding 
Vhr, the outcome transitions to the hit-and-run regime. 
In a hit-and-run event, the target retains approximately 
its original mass (although there may be some material 
exchange between the two bodies) but the escaping pro- 
jectile may be disrupted or intact. Projectile disruption 
is calculated for the reverse impact situation: the inter- 
acting mass of the target onto the total projectile mass. 
The interacting mass of the target is approximated by the 
overlapping cross sectional area of the projectile times 
the chord length through the target. The rest of the 
target mass is assumed to escape completely and is ne- 
glected in the reverse calculation. The same procedure is 
used to determine the disruption criteria for the reverse 
problem as described above; the full derivation of the re- 
verse impact is presented in LS12. When the projectile 
mass is less than about 10% of the target, the projectile 
experiences disruption during most hit-and-run events. 

For grazing collisions with impact velocities approach- 
ing V* , the collision outcome transitions from hit-and- 
run to erosion of the target. For such grazing collisions 
in the disruption regime, the outcome is defined by the 
catastrophic disruption criteria Q'^ D and Equation [13] or 
[T4l only for M h . < M t . 

In the disruption regime, recent hydrocode simulations 
of collisions between differentiated bodies showed that 
the core is preferentially incorpor ated into the largest 
remnant arcus et^[200l[20Tfl . Consider two ideal- 
ized models: 

1. Model 1 - Cores always merge: Given the orig- 
inal M corei t and M core , p , the post-impact core is 

A/corc.lr 

min(Mi r , M cor0jt + M COVCtP ). 

2. Model 2 - Cores only merge on accretion: When 



M t ). When M\ r < M t , assume that none of the 
projectile accretes and the mantle is stripped first. 
Then, A/ coro j r = min(M corCit , M\ T ). 



iMarcus et ail (|2010D found that the core mass fraction in 
the largest remnant in the hydrocode simulations fell be- 
tween these two idealized models. Hence, in this work, we 
averaged the predicted core fractions from the two mod- 
els to calculate the change in core mass fraction, / coro . 
We assumed that the initial / CO re = 1/3 to be comparable 
to Earth. 

Collision outcome maps are shown in Figure Q] for a 
mass ratio of 1:10, which is typical for the late giant 
impacts onto Earth-mass planets, and a mass ratio of 
1:40, representing the initial mass ratio of planetesimals 
colliding with embryos in the ./V-body simulations dis- 
cussed below. The impact velocity is normalized to the 
mutual surface escape velocity, and the impact angle axis 
is spaced by equa l probability according to sin 8 cos 8d8 
(|Shoemakerlll962j j. Collision outcomes are divided into 
6 groups (Figured]): 

1. Perfect merging when V < V^ sc (dark blue region 
in Figure [J) , 

2. Graze-and-merge when b > b crit and V^ sc < V < 
Vhr (blue-green region), 

3. Hit-and-run when b > 6 C rit, V; > Vh r and M\ T = Mt 
(green region), 

4. Partial accretion when b < b CI i t and M\ r > M t 
(medium blue region), 

5. Partial erosion when 0.1Af to t < M\ T < M t (e.g., 
thick black line denotes 0.5M to t), 

6. Super-catastrophic when M\ T < 0.1M tot . 

Note that collisions with impact parameters near 6 cr it 
are difficult to predict (red vertical lines in Figure [1]) , 
and outcomes near & cr it display a mix of features from 
both the accretion-to-erosion regimes and the merging- 
to-hit-and-run regimes. 

In the collision outcome maps, contours of constant fi- 
nal mass are shown by deriving the corresponding impact 
velocity as a function of impact angle, Vi i i r=cons t. Using 
the angle- dependent Q'^ D for the desired mass ratio, the 
impact energy contour is 



Q R,\r— const Q 



RD 



( M {r \ 

-2 = const 

\M tot J 



(16) 



using the universal law (Equation ITB")) for 0.1 < 
Mir/Mtot < 1. For super-catastrophic events, the im- 
pact energy contour is 



IR.\t= const 



RD 



1.8" / M lr \ 
= const 

o.i U/tct J 



l/r) 



(17) 



from Equation [TJ] Using Equation [T] the corresponding 
impact velocity as a function of impact angle is given by 



Vj lr— const 



— A/2Q_R i lr=const-Vftot/At- 



(18) 



M lT > M t , M c , 



^^COIC.t 

+ min(M core! p, Mi r 



The onset of collisional erosion is determined by the 
value of the material parameter c*. With c* = 1.9 and 
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the initial mass ratios in typical iV-body simulations, ero- 
sion begins with impact velocities of about 2 to 3 times 
the mutual escape velocity for head-on collisions and in- 
creases with impact parameter (Figure [lj . For collisions 
between similarly sized bodies, the transition between 
accretion and catastrophic disruption occurs over a very 
small range in impact velocity. For example, the critical 
impact velocities required to begin eroding and to catas- 
trophically disrupt an embryo by a body half its mass 
are only 2.0V CSC and 2.5V^ SC at b = 0, respectively. As 
the mass ratio becomes more extreme, the impact veloc- 
ities required to reach catastrophic disruption increase 
significantly. In Figure [TA, erosion begins at 2.5V CSC and 
disruption at 5.9\4 SC for a mass ratio of 1:10. For the 
1:40 mass ratio shown in Figure [TJ3, the disruption of 
an embryo by a planetesimal would require an impact 
velocity of at least 18T4 SC . 

The contours of constant largest remnant have a dif- 
ferent shape for the reverse impact in the hit-and-run 
regime compared to the forward scenario. For impact 
parameters near & cr it and M p /M t < 0.1, the interacting 
mass of the target is approximately equal to the pro- 
jectile mass. As the impact angle increases, the reverse 
projectile-to-target mass ratio decreases. Hence, the to- 
tal interacting mass in the reverse impact decreases with 
increasing impact angle. The velocity contours corre- 
spond to constant remnant mass divided by total inter- 
acting mass. The changing total mass and mass ratio 
leads to the intersection of the catastrophic disruption 
contour (blue dot-dashed lines in Figure [T]) and onset of 
projectile erosion contour (blue dotted lines) near 6 cr j t 
and divergence at larger impact angles. Futhermore, 
there is a minimum in the projectile erosion contour at 
an optimum interacting mass fraction from the target. 

3. ANALYSIS OF TV-BODY PLANET FORMATION 
SIMULATIONS 

We examined the collisions in two recent iV-body stud- 
ies of the late stag es of terrestrial planet formation. 
IQ'Brien et al.l ()2006D performed 8 simulations of terres- 
trial planet formation beginning with 25 Mars-mass em- 
bryos and an equal mass of material in a population of 
1000 .00233-MEarth p l anetes imals in a zone from 0.3 to 
4 AU. iRavmond et al.l (|2009l ) performed 40 simulations 
beginning with 85-90 planetary embryos (0.005 to 0.1 
M Earth ) and 1000 or 2000 0.0025-M Ea rth planetesimals 
between 0.5 and 4.5 AU. 

The beginning stage of each study approximates the 
end of oligarchic growth and the onset of stochastic 
growth. When the total mass in planetary embryos is 
equal to the total mass in smaller planetesimals, viscous 
stirring by the embryos cannot be damped by dynamical 
friction from the p lanetesimals, which marks the end of 
oliga rchic growth (jKokubo & Idal 119981: IGoldreich et al.l 
2004). It was assumed that the nebular gas had dissi- 
pated and Jupiter and Saturn were fully formed at the 
start of the simulations. The orbital configurations of 
Jupiter and Saturn were varied over a plausible range to 
consider their influence on the final state of planets in the 
terrestrial zone. Not all cases result in terrestrial plan- 
ets in agreement with our Solar System; moreover, the 
dynamical history of the giant planets in our Solar Sys- 
tem is still an active area of research with new ideas that 
have not yet been incorporated into detailed accretion 



studies of the te rrestrial planets (e.g. JWalsh et al.ll201 lh 
Morbidclli 120101 ) . Hence, the set of simulations are rep- 
resentative of the general dynamics of terrestrial planet 
growth in the pr esence of two outer g as giant planets. 

The study by IQ'Brien et all (f2 006) us ed the SyMBA 
jY-bo dy code ([Duncan et al.lll998f ) . an d (IRavmond et all 
120091) used the Mercury iV-body code (|Chamberslll999l) . 
Both codes are similar symplectic integrators that cal- 
culate close encounters between planetary bodies. For 
every collision, the two bodies were merged and linear 
momentum was conserved. In these simulations, the 
embryos gravitationally interacted with all other bod- 
ies. The gravitational influence of planetesimals upon 
embryos was calculated; however, planetesimals did not 
influence other planetesimals and could not collide with 
each other. In this work, we present a retrospective anal- 
ysis of the simulations to assess the range of true collision 
outcomes during the end stage of terrestrial planet for- 
mation. 

3.1. Collision outcomes 

We divided the collision outcomes from the iV-body 
studies into three groups, described in Table [TJ Group 
1 consists of all collisions leading to the growth of 15 
planets with final mass es of 0.74 to 1 . 58 Mp a rth from 8 
different simulations in IQ'Brien et alj (120061). Th ese 15 
planets were also studied bv lNimmo* et all POIOD to in- 
vestigate the evolution of the hafnium-tungsten isotopic 
system during planet growth. Group 2 includes all colli- 
sions leading to the growth of 52 planets with final masses 
between 0.70 a n d 1.4 5 Msarth from 40 simulations by 
IRavmond et al.l (|2009l ). Group 3 considers only the gi- 
an t impacts onto a l l 161 planets from the 40 simulations 
in IRavmond et al] (|2009D . Giant impacts are defined as 
collisions between planetary embryos. A final planet ex- 
periences at least one giant collision; hence, group 3 ex- 
cludes surviving embryos that only accrete planetesimals 
(26 embryos survived without a giant impact in 40 sim- 
ulations). 

In each case, the simulation collision parameters were 
used to calculate the outcome based on our analytic 
model. The outcome depends on the mass ratio of the 
two bodies, the impact angle, the impact velocity normal- 
ized to the mutual escape velocity, and the catastrophic 
disruption criteria. Here, we used values of c* = 1.9 and 
fl = 0.36 to calculate the catastrophic disruption crite- 
ria, which are appropriate for strengthless planets as we 
assumed that the planets are hot and possibly partially 
molten during the late stage of planet formation. 

3.1.1. Group 1 

There were a total of 1207 collisions by planetesimals 
and embryos during the growth of the 15 planets in group 
1. Of these, 1140 collisions were by planetesimals. The 
analytic model predicts that the majority of planetesimal 
collisions will lead to partial accretion (73%). The dom- 
inance of accretionary events is foremost a reflection of 
the impact velocity distribution, which is peaked between 
1 and 2V esc (Figure H]A-). Note that only half of the plan- 
etesimal is accreted onto an embryo when Vi = 2V esc (the 
critical velocity for accretion of half the projectile rises 
slightly with increasing impact angle, see Figure [TJ3). 
The distribution of accretion efficiencies ((Mi r — M t )/M p ) 
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Table 1 



Predicted collision outcome statistics from recent TV-body simulations of the end stage of terrestrial planet formation using the new 
collision physics model with material parameters c* = 1.9 and p, = 0.36. The TV-body simulations began with a population of planetary 
embryos and non-interacting planetcsimals and all collisions resulted in perfect merging. Giant impacts are collisions between embryos. 
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O'Brien et al. 
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Raymond et al. 2009 




Raymond et al. 2009 
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Figure 2. Distribution of impact velocities normalized to mutual surface escape velocity for (A) all collisions (black line) and all giant 
impacts (filled grey) in group 1; (B) all collisions (black line) and all giant impacts (filled grey) in group 2; (C) all giant impacts (black 
line) and last giant impact (filled grey) in group 3 TV-body simulations. 



from Asphaug 200j|) is shown in Figure The most 
common outcome of a planetesimal encounter was accre- 
tion of > 90% of the mass, although there were a signifi- 
cant number of events where less material was accreted. 
For collisions between embryos, a smaller fraction of the 
projectile is accreted (70 to 90% for the 18 partial accre- 
tion events). 

Beginning with Mars-sized embryos, there were no 
cases of erosion of the growing planet from planetesimal 
encounters. Notably, a substantial fraction of planetesi- 
mal encounters were hit-and-run events. For the initial 
1:43 mass ratio between the planetesimals and embryos, 
the critical impact parameter is about 0.78 (51°). Given 
the probability distribution of impact angles, about 40% 
of outcomes for 1 :43 mass ratio collisions are in the graz- 
ing regime, and the probability decreases as the embryos 
grow and the mass difference increases. Of all the im- 
pacts by planetesimals during the growth of Earth-mass 
planets, about 4% were graze-and-merge and 24% were 
hit-and-run. Erosion of the planetesimal occurred in 
nearly all of the hit-and-run collisions (253 out of 269), 
and catastrophic disruption of the planetesimal occured 



in about 85% of these events (Figure [TJ3). 

Next, we consider only the giant impacts in group 1. 
Giant collision outcomes are approximately evenly split 
between partial accretion, graze-and-merge, and hit-and- 
run. Only a few of the embryo projectiles in hit-and- 
run events are eroded; in these collisions, the target and 
projectile have comparable masses (7 > 0.1) and neither 
body is disrupted. 

The largest impact velocities (> 6V^ SC ), although rare, 
are high enough for embryos to catastrophically disrupt 
each other. There are notable examples of partial erosion 
(planet CJS1.4) and super-catastrophic (planet EJS1.4) 
outcomes. The last giant impact (at 222 Myr) onto a 
1.05A/Earth planet by a O.lMEarth embryo at 3.2V CSC and 
35° resulted in erosion of about 1% of the target mass 
(FigurelU'V) . The second giant impact (at 9.7 Myr) onto a 
0.21M Earth body by a 0.147M Eart h embryo at 6.2y osc and 
26° super-catastrophically disrupted the target leaving a 
largest remnant of only O.OOSMEarth- 

3.1.2. Group 2 
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Figure 3. Distribution of accretion efficiency ((Mi r — Mt)/M p ) for collisions that result in partial accretion: (A) all collisions (black line) 
and giant impacts (filled grey) in group 1; (B) all collisions (black line) and all giant impacts (filled grey) in group 2; (C) all giant impacts 
(black line) and last giant impact (filled grey) in group 3 JV-body simulations. 



There were a total of 3686 collisions by planetesimals 
and embryos during the growth of the 52 planets in group 
2. Of these, 3142 collisions were by planetesimals. As 
in group 1, the analytic model predicts that about 70% 
of planetesimal coll isions will lead t o par tial accretion. 
The simulations by iRavmond et aLI (|2009f ) considered a 
wider range of dynamical configurations for Jupiter and 
Saturn, which produced a slightly wider distribution of 
impact velocities (Figure [2)3) compared to group 1. In 
addition, the initial masses of the embryos was smaller. 
As a result, a couple of percent of planetesimals impacts 
led to erosion of the growing planet. For collisions in the 
partial accretion regime, the mean accretion efficiency is 
slightly lower in group 2 and the tail of low efficiency 
events is more pronounced than in group 1 (Figure [3)3). 

Overall, the probabilities of different collision outcomes 
for planetesimal impacts are similar in groups 1 and 2 
because of the similar mass ratios and impact velocity 
distributions (Table [1]). Again, most of the planetesimal 
hit-and-run events result in catastrophic disruption of 
the projectile. 

Compared to the giant impacts in group 1, group 2 
giant impacts have more partial accretion events and 
significantly more embryos are eroded in hit-and-run 
events. The difference is primarily a result of the fact 
th at the initial embryos were smaller in the simulations 
by lRaymond et al1(|2009f ). The larger mass ratio between 
the embryos and growing planet leads to more cases of 
fragmentation of the smaller body and fewer grazing im- 
pacts. 22% of giant impacts in group 2 have 7 < 0.1, but 
only 4% of group 1 giant impacts have such a large mass 
contrast (Figure |4|). 

Of the 151 hit-and-run giant impacts in group 2, 75 
projectiles were eroded (50%), and 29 projectiles suffered 
catastrophic disruption level fragmentation (19%). The 
three erosive giant impacts in this group removed 16, 6, 
and 1% of the material from target bodies with initial 
masses of 0.08, 0.50 and 0.72ME ar th, respectively. The 
erosive events all occurred in the eccentric Jupiter and 
Saturn (EJS) group of simulations with impact velocities 
between 2.1 and 3.2V cac . 

Because partial accretion is the most common outcome 
of non-grazing collisions, a significant fraction of giant 
impacts result in potentially observable changes in the 
bulk composition of a planet (Table [1]). For collisions 
between differentiated bodies, the core-to-mantle mass 



ratio changes during both partial accretion and erosion 
events. For the large final planets in group 2, a 5% or 
greater increase in the mass fraction of the core, / CO rc, 
occured in about 41% of all giant impacts and in about 
20% of last giant impacts. 

3.1.3. Group 3 

Considering all 1165 giant impacts onto 161 planets 
in 40 planet formation simulations, the outcomes are ap- 
proximately evenly split between partial accretion, graze- 
and-merge, and hit-and-run. Erosive events, including 
super-catastrophic disruption, occur about 1% of the 
time. In the hit-and-run events, about half of the pro- 
jectiles are eroded and about 20% are catastrophically 
disrupted. 

The equal likelihood of partial accretion, graze-and- 
merge, and hit-and-run is a result of the range of mass 
ratios and impact velocities for giant impacts. Given the 
impact velocity distribution in group 3, the velocity axis 
of a collision outcome map may be scaled by probability. 
In Figure [5[ both axes are scaled by probability; hence, 
the area of each collision outcome (denoted by colors) is 
directly proportional to their probability. 

The 4 panels span the range of giant impact mass ra- 
tios in group 3 (99% of events have 7 > 0.03 but a few 
embryo-embryo collisions have mass ratios as extreme as 
1:55). Note that, over the course of the entire simulation, 
the mass ratio of a giant impact is about equally likely 
to fall anywhere between 0.03 and 1 (Figure EC). Even 
when considering just the largest final planets (group 2), 
giant impacts may be any mass ratio (e.g., a collision 
between two 0.4MEarth bodies). However, the last giant 
impact onto target bodies greater than about 0.8ME a rth 
are dominated by mass ratios less than 0.1 (Figure lA.2F ). 

The collision maps fully scaled by probability empha- 
size the importance of the graze-and-merge regime even 
though it is a narrow regime in absolute impact veloc- 
ity. The scaled figures also emphasize that partial accre- 
tion of the projectile is the most common outcome for 
non-grazing collisions; recall that the accretion efficiency 
peaks at about 80% for all giant impacts (Figure[3p). As 
shown in Figure El hit-and-run events occur about 1/3 of 
the time and the projectile is eroded when the projectile 
mass is less than about 10% of the target. 

Note that the boundary between graze-and-merge and 
the adjacent partial accretion and hit-and-run regimes 
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Figure 4. Distribution of projectile-to-target mass ratio for all giant 
1; (B) group 2; (C) group 3 iV-body simulations. 

derived bv lKokubo fc Gend a (2010) is generally in good 
agreement with other simulations with 7 = 1. However, 
the simulations of rubbl e pile collisions with 7 = 0.5 
using the pkdgrav code (jLeinhardt et al.l [20101) found a 
much narrower graze-and-merge regime compared to the 
boundary derived from SPH simulations of fluid bodies. 
More work is needed to understand the boundaries of the 
graze-and-merge regime and its dependence on material 
properties. 

As shown in Table [TJ the last giant impact was more 
likely to be erosive (5%) or to be a hit-and-run (37%) 
compared to all giant impacts. The dynamical stirring 
by the last remaining planets leads to higher impact ve- 
locities near the end of planet formation compared to the 
total time average. In Figure note that the distri- 
bution of impact velocites is more weighted to values 
> 1.25T4 SC for the last impact (grey filled histogram) 
compared to all giant impacts (black histogram). For 
the same reason, more of the projectiles in the last gi- 
ant hit-and-run event are eroded compared to all giant 
impacts (35 out of 60 events). 

3.2. Collisions between planetesimals 

Note that the collisions considered in groups 1-3 are 
only those that contributed to the formation of the fi- 
nal planets. At the beginning of the simulation, there 
should have been many collisions between planetesimals, 
but they were not modeled. The collision outcome 
map for impacts between planetesimals with a mass ra- 
tio of 1:2 is shown in Figure [51 Typical planetesimal- 
planetesimal impact velocities will be similar to their 
collision velocities onto the embryos (about the escape 
velocity from the embryo, Figure [5]). The surface es- 
cape velocity from an Earth-mass body is about 10 
times the escape velocity from a 0.0025AfEarth plan- 
etesimal. Erosion during collisions between comparable 
mass planetesimals begins at about 3.614 sc for c* = 5, 
a value appropriate for a variety of solid compositions. 
Hence, the most common planetesimal-planetesimal col- 
lision outcome is super-catastrophic destruction of the 
planetesimals. The probability-scaled collision map (Fig- 
ure [5J3) emphasizes the extremely destructive nature of 
planetesimal-planetesimal collisions in the presence of 
growing planets. 

As shown in Figure [SJ3, collisions more oblique than 
about 60° are hit-and-run. If planetesimal-planetesimal 
collisions had been modeled in the TV-body simulations, 
the probability of a hit-and-run would have been artif- 



impacts (black lines) and last giant impact (filled grey) for (A) group 

ically high at the beginning of the calculation because 
of the assumed starting distribution of equal-mass plan- 
etesimals. In fact, planetesimals in a size distribution 
defined by a collisional cascade would have a smaller frac- 
tion of mutual hit-and-run events. The ultimate fate of 
the smallest fragments in the collisional cascade is de- 
termined by the competition between accretion onto the 
growing planets and removal from the planet's feeding 
zone (e.g., via Poynting- Robertson drag). 

In the TV-body simulations, the planetesimals all had 
the same mass because of computational limitations that 
restrict the total number particles. A tractable number 
of particles (few thousand) was insufficient to resolve a 
size distribution of planetesimals. However, the fraction 
of planetesimal collisions onto embryos that lead to par- 
tial accretion vs. hit-and-run depends on the mass ratio. 
If more of the mass in planetesimals were in smaller bod- 
ies, then accretionary collisions would be more common 
than found in the -/V-body simulations considered here. 

4. A MONTE CARLO PLANET GROWTH MODEL 

The restrospective analysis of TV-body simulations pre- 
sented above provides a limited view into the role of more 
realistic collision physics on planet growth because it can- 
not assess the cumulative effects of different collision out- 
comes. Studying suites of Monte Carlo simulations of 
the growth of a single planet via giant collisions allows 
for an intermediate examination of the role of collisions 
that is more tractable than many new full TV-body sim- 
ulations. The purpose of this Monte Carlo simulation 
is to investigate the effects of the collision physics model 
and is not meant to replace full TV-body planet formation 
models or more sop histicated statistica l population syn- 
thesis m odels (e.g.. lAlibert et al.l 120111 : iMordasini et al.l 
I2009allbl : llda fc Linll2010D . 

4.1. Method 

Our planet growth model assumed that the overall dy- 
namics of g iant impacts is s imilar to the TV-body sim- 
ulations by iRavmond et al.l (|2009fl . Because planetary 
embryos become dynamically isolated when their masses 
are much smaller than an Earth-mass planet (roughly 
one tenth the mass), most of the mass in large terres- 
trial planets, here defined a s a final mass > 0.7M p, a rth, is 
accreted via giant impacts (|Kokubo fe; Idal[l998f) . Thus, 
the Monte Carlo simulation uses the distribution of im- 
pact parameters from the group 3 TV-body simulations 
to model planet growth (summarized in Figures lA.il and 
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Figure 5. Collision outcome maps with area scaled by outcome probability. Collision maps for planetary embryos, using material 
parameters c* = 1.9, p, = 0.36, for projectile-to-target mass ratios of (A) M p : Aft = 1:1, (B) 1:2, (C) 1:10, and (D) 1:30. Vertical 
axis is scaled by velocity distribution of giant impacts in group 3 ./V-body simulations (Figure [2}. Colors are the same as in Figure [T] 
Dotted black curve - 50% of projectile accreted; dashed black curve - 80% of projectile accreted; dot-dashed blue curve - catastrophic 
disruption of projectile in a hit-and-run; dashed blue curve - onset of erosion of projectile in a hit-and-run. Example results from numerical 
simulations: fill e d circ les are merging, open circles and graze-and-merge, diamons are hit-and-run. Black - plan etesimal collisions from 
ILeinha rdt ct al. (2000); blue - Haumea formati on simulatio ns and favored Haumea impact scenario (triangle) from Lcinhardt ct al. (2010); 
arrow - disk-origin scenario for Charon from Camrp (2005); square - typical moon formation scenario from Canup (20Q4|)^ 



[Oil . In the group 3 planets, the fraction of mass ac- 
creted via planetesimal collisions varied from to 48%. 
The distribution of mass accreted from planetesimals is 
sensitive to the final mass of the planet; the large plan- 
ets accreted between 2 and 27% of their hnal mass from 
planetesimals in the stochastic stage of planet formation 
modeled in the A-body simulations (Figure lA.lj) . 
For each planet, the following procedure was applied: 

1. Randomly choose a value from the probability dis- 
tribution of initial masses (Figure [A. 1|) . The mini- 
mum initial embryo mass was O.OlMEarth- 

2. Randomly choose a value from the probabilit y dis - 
tribution of number of giant impacts (Figure [A. 

3. For each giant impact, randomly choose the projec- 
tile mass, impact velocity, and impact angle. Pro- 
jectile mass and impact velocity distributions are 
dependent on target mass (Figure 

4. Three sets of simulations were run: (i) perfect 
merging for all collisions for comparison to the 
A-body results; (ii) collision physics assuming no 
re-impact by the projectile in hit-and-run events 
(group A); and (iii) collision physics assuming 



re-impact by the projectile in hit-and-run events 
(group B). 

5. Add a randomly chosen value from the probability 
distribution of mass contributions from planetesi- 
mals. The planetesimal contribution distribution 
is also dependent on target mass (Figure IA. 1[) . 

For each simulation set, we modeled the growth of 
200 planets. When using the collision physics model, 
we tracked the core mass fraction assuming that all em- 
bryos have an initial value of 1/3 to be comparable to 
the bulk Earth. The core mass fraction of planetesimals 
was assumed to be the same as the initial embryos. We 
recorded the mass of debris produced by each giant im- 
pact. For the purpose of comparing perfect merging and 
the collision model, we assumed that all debris was lost 
and not reaccreted later. This end member assumption 
represents the maximum effect of fragmentation on in- 
hibiting planet growth in the giant impact stage. 

The distributions of impact velocities and projcctilc- 
to-target mass ratios are significantly dependent on the 
target mass (see Figure IA.2I) . The population of initial 
embryos collides with each other to grow the first mid- 
size planets. The most common collision mass ratio de- 



10 Stewart & Leinhardt: Diverse Giant Impact Outcomes 



Impact Angle [deg] Impact Angle [deg] 

30 45 60 90 30 45 60 90 




Impact Parameter b Impact Parameter b 



Figure 6. Collision outcome maps for collisions between planetesimals with projectile-to-target mass ratio M p : Mt = 1:2 and material 
parameters of c* = 5 and p, = 0.37, which are appropriate for a wide range of solid compositions. (A) Map with linear impact velocity 
normalized to planetesimal escape velocity and impact angle scaled by probability, and (B) map with area scaled to outcome probability 
for typical collision velocities stirred up by the presence of growing planets. Colors are the same as in Figure [l] Thick black curve — 
critical velocity for catastrophic disruption (M\ T = 0.5Mt o t); dashed grey curves - 10% and 90% of target mass in largest remnant; dotted 
black curve - 50% of projectile accreted; dot-dashed blue curve - catastrophic disruption of projectile in a hit-and-run; dashed blue curve 
- onset of erosion of projectile in a hit-and-run. 



pends on the initial size distribution of embryos assumed 
in a particular TV-body simulation. During this early 
period of growth in the group 3 simulations, the mass 
ratios of the two bodies was about equally distributed 
between < 7 < 1 (Figure IA.2D ). At the same time, the 
impact veloc ities are strongly peaked just above 1V CSC 
(Fi gure fA.21 A) and few collisions lead to erosion (Fig- 
ure IA.2G V As the planets continue to grow, impacts 
onto the larger bodies are dominated by smaller embryos, 
and the distribution of mass ratios is strongly peaked 
with 7 < 0.2 for Earth-mass targets (Figure IA3F ). The 
largest planets are strong perturbers on the remaining 
small embryos and the distribution of impact velocities 
becomes wider compared to earlier planet growth, span- 
ning 1 - 3V esc (Figure E2P). 

If these mass-dependent effects were not included and 
the mean distributions were utilized instead, then the 
Monte Carlo model would predict too many very large 
planets (> 1.5A/Earth)- In other words, the very end 
of terrestrial planet growth has distinct dynamical dif- 
ferences from earlier parts of the stochastic giant im- 
pact phase. To account for the variations with time, the 
values for impact velocity and mass ratio were chosen 
from the subset of group 3 collisions within ±0.1AfE a rth; 
when the target had a mass greater than one Earth mass, 
the distribution was based on the group 3 data with 
M t > 0.9M Earth . 

The projectile mass may be eroded during hit-and-run 
events. Therefore, in the group B set, we calculated the 
mass of the largest projectile remnant using the reverse 
collision scenario. The returning projectile had the mass 
of the largest remnant and the rest of the projectile de- 
bris was neglected. The re-impact had a randomly cho- 
sen impact angle and an imp act velocity given by th e 
greater of V% sin 9 or V csc after IKokubo fc tiendal l|2010h . 
The monotonically lowering re-impact velocity is an ide- 
alization that assumes that no interactions with other 



embryos led to an increase in the re-impact velocity. A 
projectile may hit-and-run several times before finally 
accreting or being disrupted. The hit-and-run sequence 
ended when the collision outcome was partial erosion, 
partial accretion, merging, or graze-and-merge. If the 
hit-and-run sequence eroded the projectile to less than 
one third of its original mass, it was assumed to be debris 
and neglected. As material was stripped from the escap- 
ing projectile, it was assumed to be derived from the 
projectile's mantle, which raised its core mass fraction. 
In this way, the core mass fraction of a planet may be 
enriched after a sequence of hit-and-run events followed 
by a merging event. 

4.2. Results 

The distribution of final planet masses was calculated 
for each simulation set. The random number seed was 
the same for each group, so the differences are entirely 
a results of the collision model assumptions. The Monte 
Carlo simulations results are summarized in Table [2] and 
Figure [3 Details of the sim ulation par ame ters i n each set 
is given in appendix figures IA.31 IA.41 and IA.5I Here, we 
focus on the properties of the largest planets with final 
masses > 0.7 A/Earth- The time between giant impacts 
was not considered here because the collision physics 
model may change the time scale for planet growth (see 

In Figure [7J the mass distribution of final planets in 
the perfect merging simulation is similar to the group 3 
A-body results, although the Monte Carlo model does 
not produce the same number of mid-size planets. The 
distribution of final planet masses is significantly smaller 
when the collision physics model is included and hit-and- 
run returns are neglected (group A). In the group 3 set 
of ./V-body simulations, 52 of the 161 planets (32%) had 
final masses > 0.7AfEaxth- In the perfect merging Monte 
Carlo simulation, 50 of 200 planets (25%) reach final 
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Table 2 

Monte Carlo planet growth simulation results for 200 planets using the collision physics model and assuming (group A) no re-impact by a 
hit-and-run projectile and (group B) with re-impact. Assuming perfect merging only, 1455 giant impacts grew 200 planets with final 
masses up to 2.0Af Bart i 1 . Of these, 50 planets with masses > 0.7M Earth grew from 651 giant impacts. For comparison, in the group 3 
TV-body simulations, 52 of 161 planets had final masses > 0.7TV/ Earth . Collision material parameters were c* = 1.9 and p, = 0.36 for 
hydrodynamic planets. 
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masses > 0.7 MEarthj however, in group A, the number 
of large planets drops to only 15 (7.5%). When hit-and- 
run return impacts are considered (group B), the final 
mass distribution of planets is between perfect merging 
and group A. In this case, 34 large planets are produced 
(17%). 

During collisional growth and fragmentation, material 
is preferentially lost from the silicate mantle, thus rais- 
ing the core mass fraction (Figure 03, E). In simulation 
groups A and B, the maximum core mass fractions are 
0.87 and 0.96, respectively, in bodies that experienced 
catastrophic impact events. Such core-dominated bodies 
are rare, and most (90-95%) final core mass fractions fall 
in the range of 0.33 to 0.4. In other words, most of the 
iron enrichment is within 20% of the initial value of / CO rc- 
However, as a group, the largest planets are more likely 
to be enriched in core mass fraction compared to smaller 
planets. In both group A and B simulations, about 2/3 
of the largest planets have core mass fractions greater 
than 10% of the initial value, compared to about 1/3 of 
all planets (Table [5]) . The largest planets experience a 
larger number of collisions which results in more cumu- 
lative erosion of the mantle. 

The mass of debris produced during planetary growth 
by giant impacts can be significant. While there were 
planets that suffered only merging collisions that pro- 
duced negligible debris (34 planets in group A and 25 
in group B), they all had final masses of less than 
0.38-MEarth and an average of only 2 giant impacts. For 
comparison, the mean number of giant impacts was 7 and 
9 for all planets in groups A and B, respectively. During 
the growth of large planets, debris production averaged 
11% of the final planet mass in group A and 15% in group 
B (Figure[7p,G). The mass of debris reported in Figure|7] 
only includes debris from giant impacts; planetesimal col- 
lisions would also have contributed to the debris during 
planet growth. 

In group A, the growth sequence that produced the 
most debris (0.28ME ar th) suffered a penultimate ero- 



sive giant impact on a planet with final mass of only 
0.22AfEarth- The most debris produced from one of the 
largest planets was 18% of a 1.31MEarth planet. No- 
tably, there is a case of a 1.2AfEarth planet in group B 
that produced 0.83ATEarth of debris during its growth. In 
some cases, the growth sequence includes a step where 
the largest remnant is smaller than the initial embryo. 
Such destructive sequences occured for 4 planets in group 
A and 7 planets in group B. 

An average of 12 giant impacts grew the largest plan- 
ets in group A. With the inclusion of hit-and-run re- 
turn events, an average of 16 giant impacts grew the 
largest planets in group B. In group B, the largest num- 
ber of giant impacts for any planets was 26 (an ultimately 
0.89MEarth planet), in contrast to the maximum of 18 
giant impacts in group A. Hit-and-run return collisions 
often led to multiple re-impact events before the final 
merging or disruption of the projectile. The number of 
excess giant impacts in group B is shown in Figure El 
The mean number of extra collisions on the largest plan- 
ets was 4 and ranged from to 8. 

5. DISCUSSION 

This work demonstrates that the final, stochastic stage 
of terrestrial planet formation encompasses a diversity of 
collision outcomes. All types of collisions, from super- 
catastrophic disruption to perfect merging, are possible. 
Previous work that assumed perfect merging for all colli- 
sions were not capturing the full complexity of the giant 
impact phase. For the expected dynamical conditions 
for the terrestrial planets in our early Solar System, the 
principal collisions outcomes are approximately evenly 
split between partial accretion, hit-and-run, and graze- 
and-merge. 

Because the collision outcome depends on mass ratio, 
the initial growth of planets via giant impacts is sensitive 
to the assumed mass distribution of planetary embryos. 
Because different numerical methods are typically used 
to model the growth of embryos (the oligarchic phase) 



12 



Stewart & Leinhardt: Diverse Giant Impact Outcomes 



Oh 



g 

Ph 




0.0 0.5 1.0 1.5 
Final Mass [M Earth ] 



| 
o- 



as 




0.0 0.5 1.0 
Final Mass [M Earth 

0.4 




I 
a, 



CO 
Ph 



M 

0- 

o 

c 

o 



p- 




0.0 0.5 1.0 1.5 
Final Mass [M Earth ] 



M 

Oh 

O 

C 

O 



P- 




0.5 

0.4 

0.3 

0.2 

0.1 
0.0 



M debris /M final 



0.0 0.2 0.4 0.6 0.8 1.0 

Mdebns /M fma , 



Figure 7. Monte Carlo planet growth simulation results. Mass distributions, core mass fractions, and debris mass (grey filled histograms) 
for 200 planet growth calculations assuming (A) perfect merging, (B,D,F) collision physics model with no hit-and-run return collisions, and 
(C,E,G) collision physics model with re-impact by hit-and-run projectile. Black line histogram in A-C is mass distribution of planets from 
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and the end stage of planet formation (the stochastic 
phase), the discontinuity in technique has led to sig- 
nificant variation in the assumed initial number of em- 
bryos spread over different radi al distances in d i fferent 
iV-body studies. For example, iRavmond et al.l (2009) 
beg an with 85 to 9 em bryos between 0.5 and 4.5 AU 
and lO'Brien et al.1 (|2006[ ) assumed 25 embryos between 
0.3 and 4 AU. Both include d a population of equ al-sized 
planetesimals. In contrast, iKokubo et al.l (|2006f ) began 
with 5 to 30 embr yos between 0.5 and 1.5 AU and no 
planetesimals, and feokubo fc Gendal (|2010[) assumed 16 
initial embryos with no planetesimals over the same dis- 
tance from the Sun. 

In all ./V-body simulations, the initial masses and ra- 
dial distributions of embryos are motivated by the isola- 
tion mass achieved a t the end of oligarchic growth (e.g., 
IKokubo fc Idalfl998l) . However, the time scale to reach 
isolation is a function of radial distance and surface den- 
sity, whereas the initial state of iV-body simulations im- 
plies that the isolation masses have been reached every- 
where at the same time. As a result, the number and 
size distributions of embryos and planetesimals are not 



perfectly linked to earlier stages of planet formation. 

To address this problem, a few groups have recently de- 
veloped hybrid computational techniques that may link 
more seamlessly the oligarchic gro wth and the stochas- 
tic end stage of planet format i on dLevison et al. 120051 : 
iBromlev fc Kenvonl [2001 120 lit iMorishima et aljl2010f ). 
The current suite of hybrid techniques are based on very 
different approaches with different strengths and weak- 
nesses. In this work, we have demonstrated that the end 
stage of planet formation involves significant fragmenta- 
tion. Because the number of gravitationally interacting 
particles is limited in direct iV-body techniques, hybrid 
methods that may include large numbers of smaller bod- 
ies (dust to planetesimals) are required to be able to fully 
model the giant impact stage. 

5.1. Fragmentation during the giant impact stage 

Fragmentation is a significant component of the giant 
impact stage of terrestrial planet formation. Fragmenta- 
tion occurs primarily during partial accretion events and 
in erosion of the projectile in hit-and-run events. About 
1 /3 of giant impacts lead to partial accretion and about 
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1/3 are hit-and-run events. Of the hit-and-run events, 
about half are energetic enough to erode the projectile 
(Table [2]). A small amount of escaping debris is created 
in some graze-and-merge events; however, the amount 
has not been quantified in our simulations as it is com- 
parable to the resolution lim it (of order 1% of the mass) 
i|Leinhardt fc Stewart! I2012T) . In Figure the contours 
for projectile erosion (blue dashed line) extend into the 
graze-and-merge regime for collisions with a mass ratio 
more extreme than 1:7. Overall, about half of all giant 
impacts create significant (greater than about 1% of the 
total mass) debris. 

In contrast, events that lead to erosion of the larger 
body are rare. All types of target erosion events (partial 
erosion to super-catastrophic disruption) occur about 1% 
of the time. Although infrequent, the conditions required 
to strip Mercury's mantle by a single giant impact are 
included in the range of collision outcomes in the early 
Solar System. 

To date, numerical simulations that include full frag- 
mentation models have focused on earlier stages of 
planet forma tion, from dus t to planetesimals to oli- 
garch s (e.g.. iKenvon fc Luul 119991: iBromlev fc Kenvon 
2011 1 : ILeinhardt fe Richardson! 120051: ILeinhardt et al. 
20091 iChambers! 120081: iKobavashi et al.N2010ft . In these 
studies, fragmentation laws were defined by a size- 
dependen t catastrophic dis r uptio n criteria such as de- 
rived by iBenz fc Asphaugl (|1999f ) (similar to equation 
[21 but including a strength regime for bodies smaller 
than about 1 km in radius). The previous fragmenta- 
tion models have assumed that the disruption criteria 
follows pure energy scaling. Under pure energy scaling, 
the magnitude of the impact velocity does not influence 
the threshold energy for disruption (and \x — 2/3 in equa- 
tion^. In other words, the disruption criteria for a par- 
ticular size target only dep ended on the specific energ; 



ic energy 
~M C 2 009) 



of t he impact M n V„ ■ /M t . InlStewart fe Leinhari .. 
and ILeinhardt fe Steward (|2012f ). we showed that this as- 
sumption is invalid. 

In fact, catastrophic disruption f ollows nearly pure mo- 
ment um scaling, where jl — 1/3 (|Leinhardt fe Stewart! 
2012). As a result, the disruption threshold is sensitive 
to the magnitude of the impact velocity by the factor 
y*{2-3p.) j n equation[2] Because most numerical and lab- 
oratory disruption experiments have involved projectiles 
that are much smaller than the target, most of the pub- 
lished disruption criteria are specifically for high impact 



velocities and extreme projectile-to-target mass ratios. 
Hence, when the collision involves two bodies that are 
more similar in mass, the applied disruption criteria are 
too high and the frequency of disruption is underesti- 
mated; 

IStewart fe Leinhardt! ()2009[ ) compiled disruption data 
over a wide range of mass ratios in both the strength 
and gravity regimes to conclude that energy scaling was 
incorre ct and the results were close r to momentum scal- 
ing. In ILeinhardt fc Stewart! (|2012l ). we derived a general 
analytic model for the disruption criteria in the gravity 
regime and fit material parameters for a range of tar- 
get body types from solid planetesimals to fluid planets. 
Future work will tackle catastrophic disruption in the 
strength regime. 

Now, future models of planet formation have a robust 
analytic model for fragmentation during collisions be- 
tween gravity-dominated bodies. Our estimate of the 
magnitude of debris production, on average 15% of the 
mass in the largest planets, is in excellent agreement 
with preliminary results from multi-scale A -body and 
hydrocode calculations. iGenda et al.l (|2011l ) calculated 
the very end stage of planet formation, beginning with 
16 embryos and a total mass of 2.3 MEarth, using an 
A-body code. For each collision, the outcome was calcu- 
lated by an SPH hydrocode simulation. The largest rem- 
nant or the two hit-and-run bodies were re-inserted into 
the TV-body calculation and the debris was neglected. 
The cumulative mass of the fragments was about 22% 
of the total mass (0.48MEarth)- Such multi-scale calcula- 
tions are extremely computationally expensive, and the 
analytic model will allow for more and more detailed in- 
vestigations of planet formation. 

The need to track many small fragments is a signifi- 
cant challenge in models of the end stage of planet for- 
mation. The dynamical interactions between the proto- 
planets must be calculated by direct TV-body techniques, 
but the smaller fragments will need to be treated statis- 
tically. As a result, detailed models of the giant impact 
stage must adopt new methods, such as the hybrid codes 
mentioned above. Hybrid models are needed to be able 
to calculate what fraction of the debris is reaccreted onto 
the final planets and what fraction is removed (e.g., via 
Poynting- Robertson drag). For fragments ground down 
below about 1 km in size, new disruption criteria are still 
needed in the strength regime that fully account for ma- 
terial properties, impact angle, mass ratio, and impact 
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velocity. 

5.2. Time scale of planet formation 

The influence of more realistic collision physics on the 
time scale of terrestrial planet formation is difficult to 
predict because of competing factors. On one hand, the 
growth rate of planets is slower when outcomes other 
than perfect merging are included. On the other hand, 
the debris produced by collisions can influence the overall 
dynamics of the embryos and planetesimals (e.g., via dy- 
namical friction, which could lead to more lower velocity 
collisions and more frequent merging outcomes). 

The time scale for the end stage of planet formation 
was investigated in a recent study that in cluded two col- 
lision outcomes. iKokubo fc Gendal (|2010f ) conducted N- 
body simulations with a collision model that allowed for 
either perfect merging or an ideal hit-and-run event us- 
ing their empirical boundary presented in equation [15] 
(which they applied at all impact angles). In a hit-and- 
run collision, neither body lost mass (there was no frag- 
mentation) but the relative velocities of the bodies de- 
creased. In simulations that began with 16 equal-mass 
embryos between 0.5 and 1.5 AU, they found that about 
half of the collisions were hit-and-run events. After a 
hit-and-run encounter, the reduced relative velocities led 
to a high probability of merging on the subsequent en- 
counter. As a result, the time scale for planet growth 
was essentially the same as in simulations with the same 
starting conditions that assumed perfect merging. 

The magnitude of any changes to planet growth time 
scales from different collision outcomes is sensitive to the 
initial conditions in the simulation. The simulations by 
IKokubo fc Gendal (|2010ft began with few relatively large 
embryos (0.15M®) and no planetesimals. With this ini- 
tial distribution of embryos, all collisions were between 
comparable mass bodies (M p /M t > 0.1), which are the 
most likely to be graze-and-merge or hit-and-run events 
(Figure [5]). Since the hit-and-run bodies are likely to 
recollide on time scales comparable to the orbital period, 
the overall effect of hit-and-run events on this very end 
stage of planet growth was found to be negligible. If 
fragmentation were included or if the sizes of the bodies 
were more diverse (e.g., with the inclusion of planetesi- 
mals or smaller initial embryos), more realistic collision 
outcomes would have a larger effect on planet growth. 
As mentioned above, using th e same initial condit ions in 
their multi-scale calculations, iGenda et al.l (|2011ft found 
that fragmentation was a significant process. 

Although fragmentation makes planet growth from 
an individual collision less efficient, other effects from 
the debris may lead to faster planet growth overall or 
faster stages of planet growth. Fo r example, during oli- 
garchic growth. IChambersl ()2008l ) found that fragmen- 
tation led to faster growth of embryos because smaller 
fragments are more easily captured. However, fragmen- 
tation also decreased the surface density of solids in the 
disk (which limits the embryo's final mass) because frag- 
ments were lost more quickly by drag processes as they 
were ground down in size . During the giant impact 
stage. lO'Brien et"aT1 ((2006) found that the strong dy- 
namical friction from 1000 planetesimals led to overall 
faster planet growth compared to studies without any 
planetesimals. However, their study assumed perfect 
merging for all collision outcomes and non-interacting 



planetesimals. If dynamical friction becomes very large, 
a gap could form in the planetesimal disk around an em- 
bryo, which would effectively halt the growth of that em- 
bryo. Hence, the role of fragmentation on planet growth 
time scales is not independent of other processes acting 
at the same time. Ultimately, because fragmentation is 
a critical process that feeds back into the dynamics of 
planet growth, new simulations that include both em- 
bryos and a fully interacting population of small bodies 
are needed to investigate how more realistic collision out- 
comes affect formation time scales. 

5.3. Composition and chemistry of planets 

We found that the core-to-mantle mass fraction in- 
creases during the growth of planets via fragmentation 
during collisions between differentiated embryos. Loss 
of the silicate mantle primarily occurred during partial 
accretion and erosion of the projectile in hit-and-run 
events. From our group B Monte Carlo calculations, 
the mean increase in the core mass fraction was 15% 
for the 34 largest planets (Figure [35, black histogram). 
The magnitude of the core fraction increase is potentially 
observable in the study of the chemical composition of 
p lanets and early Solar Sy stem materials. 

lO'Neill fc Palmd (|2008l i argue that the Earth's bulk 
iron to magnesium (Fe/Mg) ratio is significantly larger 
than the solar ratio. They propose that the Earth lost 
a portion of its silicate mantle during the giant impact 
stage, which raised the Fe/Mg ratio of the final planet 
compared to the more primitive (closer to nebular com- 
position) materia l s that formed the planetary embryos. 
lO'Neill fc Palmd (|200gl estimate that the whole-Earth 
Fe/Mg mass ratio is 2.1 ±0.1. The Fe/Mg value for prim- 
itive materials is not know n precisely. The so lar pho- 
tosphere value, 1.87 ± 0.4 (|Asplund et al.l 120051 ). is too 
poorly constrained to be useful for such detailed com- 
parisons. The Fe/Mg ratio for the solar wind will be 
constrained by Genesis mission. Early results suggested 
a lower value than the solar photosphere (1.61 ± 0.23, 
Uurewicz et al.ll2011h : however, final data calibration is 
still in progress (A. Jurewicz, pers. comm.). The Fe/Mg 
ratio for carbonaceous chondrite s, the most primitiv e 
type of meteorite, is 1.92 ± 0.08 (jPalme fc Jonesl 12005ft . 
The available data suggest that the Earth is enriched in 
Fe/Mg compared to solar composition by approximately 
10%. 

Our Monte Carlo calculations of the cumulative ef- 
fects of realistic collision outcomes are in excellent 
agreement with the coll isional erosion idea proposed by 
lO'Neill & Palmd ([2008). Although our planet growth 
models are not detailed enough to track individual el- 
ements, we use the core-to-mantle mass ratio as a proxy 
for the Fe/Mg ratio. About half of the largest planets 
in the group B calculations have core mass fraction en- 
richments of 10-30%. The mean is likely to be slightly 
lower than this range because some of the mantle frag- 
ments will be reaccreted. Assuming that the reaccretion 
process is less than 100% efficient, the largest terrestrial 
planets have a larger core mass fraction compared to the 
initial embryo composition as a result of fragmentation 
during giant impacts. 

Bulk elemental ratios are difficult to derive for a planet, 
and more detailed information can be derived from ex- 
amination of isotopic systems. The hafnium-tungsten 
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(Hf-W) system, with a half life of 9 Myr, is a major 
constraint on the timing of planet formation. Both ele- 
ments are refractory (high condensation temperatures); 
Hf is lithophile and is retained in the silicate mantle, 
while W is moderately siderophile and prefers the iron 
core. With assumptions about the extent of chemical 
equilibration between the metals and silicates within the 
growing planet, measurements of the present-day ratios 
of Hf and W isotopes constrain the time scale for planet 
formation, inclu ding the timing of the proposed M oon- 
forming impact (|Jacobsenl l2005: Kl eine et alll2009l ). 

Using 15 p l anets from the iV-body si mulations by 
lO'Brien et all (|2006l) (Group 1 in Table [Tjl. lNimmo et all 
( 20101) calculated the evolution of the Hf-W system dur- 
ing planet growth in the giant impact stage and com- 
pared the results to observed values for Earth, the Moon, 
and Mars. The study varied the equilibration factor 
between two idealized end members: no equilibration, 
where the core from the projectile merges with the target 
core without any equilibration with the mantle, to per- 
fect equilibration, where the projectile core completely 
equilibrates with the mantle. The equilibration factor 
was held constant during the entire growth of a planet; in 
reality, the extent of equilibration depends on the physics 
of material mixing during gia nt impacts, which is sti ll 
rather poorly unders tood (e.g.. lDahl fc Stevenson|[2010l ) . 
iNimmo et all (|2010f ) found that Earth-like Hf-W ratios 
could be generated if the iron cores partially equilibrated 
with the mantle; however, they could not simultaneously 
match both the Earth and the Moon. They suggest that 
either that th e Earth and Moon equilibra ted after the 
giant impact ( Pa hlevan fc Stevenson! l20"07j) or that the 
time scales for planet growth are too short in the TV- 
body simulations. 

With more realistic collision physics in planet forma- 
tion simulations, calculations of the evolution of the Hf- 
W system will be more robust. In particular, the colli- 
sion model estimates the fraction of core and mantle in- 
corporated during partial accretion of the projectile and 
following hit-and-run events with projectile erosion. Im- 
proving the time scale for planet growth, as discussed 
above, is essential for interpretation of the Hf-W system. 
Finally, improved physical models for mixing of metals 
and silicates during giant impacts are still needed. 

Fragmentation during planet formation may also affect 
the final volatile content. Several studies have investi- 
gated possible sources of water in the Earth by assum- 
ing that the water content in condensed material varies 
monotonically with distance from the Sun. The initial lo- 
cations of material that accretes into terrestrial planets 
have been used to investigate the formation of Earth-like 
planets under th e ass umption of perfe ct m erging (e.g., 
iMorbidelli etalll2000l IRavmond et alll2004 12007). The 
accretion of water to the Earth may now be investigated 
in greater detail. 

The incorporation of volatiles in planets not only de- 
pends on the source location of the in coming materia l 
but also the specific collision scenario. lAsnhauel ll20Toh 
proposed that the volatile content of planetary bodies 
may be affected by hit-and-run events. Because erosion 
of the projectile is common in hit-and-run events, the 
projectile may be stripped of volatile-rich outer layers 
(including an atmosphere). If the hit-and-run projectile 
is later incorporated into a growing planet, the planet 



would be depleted in both volatiles and mantle material 
compared to the initial embryos. 

The process of compositional changes via fragmenta- 
tion is not restricted to the inner Solar System. The 
dwarf planets in the outer Solar System have higher 
bulk densities than smaller Kuiper Belt objects (|Brownl 
l2008HFraser fc Brownll2010D . The bulk densities of dwarf 
planets may have increased during collisional growth by 
preferential stripping of the icy mantles from projectiles, 
possibly by partial accretion of planetesimals during run- 
away growth in addition to the limited number of giant 
impacts experienced by dwarf planets. 

5.4. The graze- and-merge regime 

The boundary between the graze-and-merge and hit- 
and-run regimes is not precisely known. In this 
work, we used the boundary defined from a single 
set of calculations for differentiated iron-silicat e plan- 
ets (jKoSboXGendi HSU [GendaeTaD Hoil)- The 
boundary is slightly different in the few other studies that 
have focused on this regime. In Figure[5j the results from 
diff erent studies are plotted a s symbols for comparison to 
t helKokubo fc Gendal 1 2010) graze-and-merge boundary. 

iLeinhardt et a l. (2000) calculated the outcome of colli- 
sions between equal- mass rubble piles (black symbolf0). 
The collision velocities were subsonic between the 1- 
km radius bodies, so the study utilized the pkdgrav TV- 
body code rather than a shock hydrocode. The results 
show that grazing outcomes extend to impact parameters 
slightly below & cr ;t. For impact angles between about 30 
and 50 degrees, t he transition to hit-and -run is in good 
agreement with (|Kokubo fc Gendal 12010). However, at 
higher impact angles (64°), the rubble pile collisions tran- 
sitioned to hit-and-run at much lower impact velocities. 
Using the same code, ILeinhardt et all (|2010l ) modeled 
collisions onto ~ 700 km radius bodies with equal mass 
and half mass projectiles (blue symbols in Figure [5j\ and 
B). The resul ts for equal mass bodies support the same 
boundary as (jKokubo fc Genddl20Tfj| ) for the 30 to 40° 
impacts. However, the simulations for a mass ratio of 
1:2 show a much faster transition to hit-and-run at an 
impact angle of 53°. 

The different boundaries between graze-and-merge and 
hit-and-run may be primarily attributed to the relatively 
low resolution for all three studies. The studies were fo- 
cused on other aspects of collisions and not to designed 
to resolve the interaction of very thin layers of material 
at high impact angles. In addition, small differences in 
dissipation of energy and momentum are likely to influ- 
ence the transition to hit-and-run at high impact angles. 
If the graze-and-merge regime is smaller than considered 
here, the collisions would increase the fraction of hit-and- 
run events. 

The graze-and-merge regime is believed to have left 
a strong mark in our Solar System. The form ation of 
Earth 's Moon (|Canupl l200l) . the Pluto system dCanupl 
120051 ). and the Haumea system (jLeinhardt et al.l 120101 ) 

These data arc derived from Table 1 in[Lemhardt ct al. (2000). 
The graze-and-merge regime is identified as when the accreting 
mass is equal to the projectile mass (e.g., cases where the graze- 
and-merge was a long duration event); some of the points plotted 
as perfect merging outcomes may have gone through a fast graze- 
and-merge event. 
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are all attributed to graze-and-merge regime events (Fig- 
ure [5]). More work is needed to understand the details of 
moon formation in the graze-and-merge regime and the 
factors that control the transition between the graze-and- 
merge and hit-and-run regimes. 

6. CONCLUSIONS 

Using a new analytic collision physics model 
(|Leinhardt fc Stewartl |2012| ), we have investigated the 
range of collision outcomes during the stochastic end 
stage of planet formation. For the dynamical conditions 
expected in our early Solar System, the outcome of giant 
impacts span all possible regimes: hit-and-run, merg- 
ing, partial accretion, partial erosion, and catastrophic 
disruption. Fragmentation during giant impacts is sig- 
nificant. During the formation of planets larger than 
0.7AfEarth, the total mass of debris is about 15% of the 
final planet mass. Fragmentation occurs primarily by 
erosion of the smaller body in partial accretion and hit- 
and-run events. Future simulations of the end stage of 
planet formation will need to utilize hybrid techniques 
that are able to track both massive planets and a large 
population of smaller bodies. Assuming that fragments 
are not completely reaccreted, growth via giant impacts 
creates final planets are that depleted in volatiles and 
mantle material compared to the initial planetary em- 
bryos. 
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Figure A.l. Distribution of parameters from group 3 Af-body simulations. Distributions for initial mass, number of giant impacts, and 
mass contribution from planetesimals used in Monte Carlo planet growth simulations. 
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ratios used in Monte Carlo planet growth simulations. 



The Astrophysical Journal 



o 
a 

Cl. 

o 

n 
o 



o 
Ph 

E 




0.0 0.5 1.0 1.5 
Final Mass [M Earth ] 



1.0 
3 0.8 

s 

P n -7 



Oh 

CO 
CO 



0.6 
0.5 













o ocfe 

o 


f&B o 


° So ° 


f » o 3> ° 




S o°° 


So o 


<?° °|o ° 


o 

o 


O o ° 




n O 

° o° 




° o 
° °0° % 











0.0 0.5 1.0 1.5 
Final Mass [M Earth ] 



CD 

B 

Cl. 

O 

a 

o 



Li- 




re 
CL 

s 
5 

o 

CD 
g 

3 

Z 



.00 0.05 0.10 0.15 
Initial Mass [M Earth ] 



20 
15 

10 

5 




oo oaoo o o o 
o o 

OOOCGSO QD OO© O 

o o o oo o o 

O OTD IB O <IBO O GD O 
O OCEQ SSDOBB O O CO 



0.0 0.5 1.0 1.5 
Final Mass [M Earth ] 



0.0 0.5 1.0 1.5 
Final Mass [M Earth ] 




CD 

M 0.08 \r 

<3 0.06 

d 
o 



10 15 20 
Number of Giant Impacts 




0.0 0.2 0.4 0.6 0.8 1.0 

M p /M t 




0.0 0.5 1.0 1.5 2.0 2.5 



Figure A. 3. Distribution of parameters for 200 Monte Carlo planet growth simulations with perfect merging. 
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Figure A. 4. Distribution of parameters for 200 Monte Carlo planet growth simulations with collision physics model and assuming no 
re-impact by a hit-and-run projectile (group A). 
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Figure A. 5. Distribution of parameters for 200 Monte Carlo planet growth simulation with collision physics model with re-impact by a 
hit-and-run projectile (group B). Note the increase in number of giant impacts compared to group A. 



